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We express the in-in functional determinant giving the one-loop effective potential for a 
scalar field propagating in a cosmological spacetime in terms of the mode functions speci¬ 
fying the vacuum of the theory and then apply adiabatic regularization to make this bare 
potential finite. In this setup, the adiabatic regularization offers a particular renormalization 
prescription that isolates the effects of the cosmic expansion. We apply our findings to deter¬ 
mine the radiative corrections to the classical inflaton potentials in scalar field inflationary 
models and also we derive an effective potential for the superhorizon curvature perturbation 
C encoding its scatterings with the subhorizon modes. Although the resulting modifications 
to the cosmological observables like nongaussianity turn out to be small, they distinctively 
appear after horizon crossing. 


I. INTRODUCTION 

Functional determinants arise in various instances in quantum field theory like in the calculations 
of the effective actions, gauge fixing Faddeev-Popov terms, semiclassical tunneling amplitudes and 
Jacobian factors (for a review see e.g. [1]). In general, they appear as results of Gaussian path 
integrals and there are different methods in evaluating them such as the heat-kernel expansion, 
zeta function regularization and Gel’fand-Yaglom theorem. 

The functional determinants also appear in quantum cosmology since Gaussian in-in/Schwinger- 
Keldysh path integrals are often encountered. For example, in certain models of preheating one 
assumes the existence of a reheating scalar that interacts with the inflaton, yet the action may 
still be quadratic in the reheating scalar yielding a Gaussian path integral in quantum theory 
(see [2, 3]). Similarly, in single scalar field models one may consider integrating out the quadratic 
inflaton fluctuations about the inflating background to determine the quantum backreaction effects 
(see e.g. [4-9]). Obviously, the functional determinants in a cosmological setting are time dependent 
and they can naturally be interpreted as effective actions. 
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Most of the time the physically relevant quantity is not the sole functional determinant of an 
operator but the ratio of the determinants of two closely related operators. In [14] a simple formula 
that gives the value of a normalized Gaussian path integral with a time dependent frequency is 
derived and we adapt that formula to the cosmological in-in functional determinants in the appendix 
(a formally similar trick has been frequently used in cosmological applications but the derivation 
of [14] is rigorous). It turns out that in a cosmological setting the result can be expressed in terms 
of the mode functions identifying the vacuum of the theory. This simple observation allows one 
to work out the Gaussian path integrals in several cosmological backgrounds with different vacua 
and utilize various approximation methods. We show that the bare Coleman-Weinberg effective 
potential [15] can be obtained by using the standard mode functions of the flat space. 

The effective potential obtained after a Gaussian integral involves a single loop momentum 
integral and it can be identified as a one-loop approximate result that is otherwise exact. Typically, 
this formal expression has the standard ultraviolet (UV) divergence which must be cured by a viable 
method. On the other hand, in (quasi) de Sitter space one also encounters peculiar infrared (IR) 
divergences in loop calculations (see e.g. [16, 17]). As discussed in [18], different regularizations 
agree in de Sitter space when the UV and IR cutoffs are chosen to be constants in physical and 
comoving spaces, respectively. 

The dimensional regularization has been used extensively to calculate the one-loop effective 
potentials in de Sitter space, see e.g. [4-9] (see also [10-13] for different approaches). Therefore, 
in this paper we prefer to utilize adiabatic regularization [19, 20] in curing the loop infinities 
of cosmological functional determinants. Recently, the adiabatic regularization has been used 
to renormalize the inflationary power spectrum revising some of the well known results [21-23]. 
We believe that in the absence of any direct observational data, it is useful to compare different 
regularization methods to have a clearer picture of the physics. 

As we will discuss below, the adiabatic regularization captures the contribution of the cosmic 
expansion on the quantum effective potential. Namely, it encodes “particle creation effects” rather 
than measuring “vacuum polarization”. We show that the terms subtracted by the adiabatic 
regularization can be interpreted as counterterms appearing in the bare action. Moreover, the 
calculation is IR safe, i.e. no IR divergences appear that require a careful treatment. 

In this paper, we also consider the coupling of the superhorizon curvature perturbation to 
the subhorizon modes of a scalar field in a quasi-de Sitter spacetime with constant deceleration. 
Since ^ is a special metric function, the coupling we are dealing with is fairly universal. By 
integrating out the Gaussian subhorizon modes in the path integral, it is possible to obtain an 



3 


effective potential for (^gh encoding its scatterings with subhorizon modes. We determine how 
this quantum effective potential affects the power spectrum and nongaussianity. As expected, 
the corrections turn out be suppressed by the factor H'^/Mp but they also appear after horizon 
crossing, which is a distinguishing feature. 

II. A SAMPLE PROBLEM 

Consider two real scalar fields cj) and x that are minimally coupled to gravity with the standard 
action 

S = -\f d"xV^[V^,/.V^</. + V^xV^X + 2P(0,x)], (1) 

where the potential is assumed to have the following form 

V (</), x) = v(4>) + ^m^x^ + (2) 

The scalars are taken to be propagating in a cosmological background with the metric 

ds^ = —dt^ + a{t)^ dx''dx^. (3) 

Suppose that one is interested in calculating the correlation functions oi cj)va. a, suitable vacuum of 
the theory using the in-in formalism. Since the action is quadratic in the x held, one may temp to 
integrate it out completely in the in-in path integral. Schematically written, the relevant part in 
the path integral involving the x held takes the form (see e.g. [25]) 

I Vx. Dx+ Dx- DP+ DP- exp |i j'* (P+x+ - H+) - ^| 'l/[x+(to)]^[x"(to)], (4) 

where H is the Hamiltonian and T denotes the vacuum wave-functionak of the x held defined 
at an initial time to- Since it is very difficult to determine the vacuum of the interacting theory, 
one usually takes to ^ (when it is possible to do so), and in that case the sole effect of 

the vacuum wave-functionals in (4) is to produce the necessary ie terms for the propagators in 
the perturbation theory [24, 25]. The standard in-in path integral measure amounts to summing 
over all doubled phase space fields x^ and P^, where the fields x^ are constrained to satisfy 
X~^(t*) = = xit*) and Vx* denote the integration over these x(^*) configurations defined at 

the fixed time t*- 

^ In an interacting theory, the vacuum wave-functional may not be written as the product of 4> and x pieces. However, 
this will be the case for us since we eventually examine the to ^ —oo limit. 
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One can carry out the quadratic integrals over momenta in a straightforward way. Defining 
a new field as 




(5) 


one can see that the scale factor a{t) dependent factor coming to the measure from the Gaussian 
integrations is exactly canceled out by the Jacobian of the transformation (5). As a result, the 
integral (4) becomes 

J Dii~ exp — 5“)] , (6) 

where 


S 



— {ON 


0 


±1 “T 


(7) 


In (6), the vacuum wave-functionals are suppressed since their only role is to prescribe the propa¬ 
gators of the theory. Note that, here (j) plays the role of an external field, which will be integrated 
eventually to give the cj) correlation functions. 

Defining 


and 


wehere 


M- = 




L 0 


' g\ct>+f 0 

Ls — 




0 -L 


0 _ 


L = 


dP 


1 


9 


- Tidf —-H — -H J- m 


(6) can be rewritten as 


J D/r* Djx exp 




( 8 ) 


(9) 


( 10 ) 


( 11 ) 


As discussed in [26], the in-in path integral is actually over the fields satisfying ^"’'(t*) = ^“(t*) and 
= /i“(t*)- These boundary conditions ensure the absence of surface terms that otherwise 
arise after integrating (7) by parts to yield (11). Moreover, they also make the operator 
essentially self-adjoint [26]. Since p. is a real field, the Gaussian integral (11) gives (^[det Li]“^/^. 
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The normalization C can be fixed from the requirement that the path integral yields unity when 
0 = 0, which evaluates (11) as 

det Li 
det Lq 

One can now use (A.20) to express this ratio of determinants to reach 

ds J d!^x [G++[ 0 +, 0 “](x;x) 5 r^ 0 +(x)^-G 7 "[ 0 +, 0 "](x;x) 5 r^ 0 "(x)^]|, (13) 

where and G~g are the matrix entries of the Green function G* of the operator (9) as defined 
in the appendix. Note that these Green functions depend on 0^ in a very nontrivial way. As a 
result, integrating out the y field gives the following bare in-in effective action 

5’e//[0’^,0“] = ^ y j [G'++[0+,0"](x;x)0+(x)2 - G7“[0+,0 “](x;x)0“(x)2] , (14) 

which must be though to correct the classical 0 action given in (1). 

For given 0^ fields, one can in principle solve for the Green function G^, which is assumed 
to be fixed uniquely once the vacuum of the theory is specified. However, it is difficult, if not 
impossible, to carry out this computation for arbitrary fields. Therefore, we consider constant 
0^ configurations to read the quantum effective potential. In that case, ^ becomes a free massive 
scalar field propagating in (3) and it is easy to calculate the corresponding Green functions. Indeed, 
quantizing the field by introducing the standard ladder operators and the mode functions as 

tJ-s = ^2^3/2 / + (15) 

where [ag, atj = d'^ik — k') and 

k 9 9(9 9 9 3 • 

+ + sg^(j)^ - - -H fik = 0, 

I 4 2 

MfeAfc - k*kt^k = i, (16) 

the Green functions are determined by the following vacuum expectation values of the operators 

G++ = i (Tifittit)), G~ = t (T(^7/i7)), (17) 

where T and T refer to the time and anti-time ordering operations, and the vacuum is defined as 
usual by ag|0>= 0. Using these in (14), one finds for constant 0^ that 

5e//[0+,0“] = y d^xa^ [Veffi^P^) - 14//(0")] , 




( 18 ) 
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FIG. 1. The perturbative one-loop graphs arising in the calculation of the effective potential (19). 

where the bare effective potential is given by^ 

1 r rl^k 

y (^/o (19) 

Recall that the mode functions n[t,s,(p] obey (16) and they are uniquely fixed when the vacuum 
is specified. Using (5), one can also express using the y mode functions by Hk = 

so that the initial scale dependent factor in (19) drops out. As discussed in [27, 28], the classical 
configurations arise as stationary phases of the in-in path integrals. Consequently, after a suit¬ 
able renormalization Vjjj:{(f>) must be thought to modify the tree-level potential u((/>) given in (2) 
governing the classical evolution of the (/> field. 

The effective potential Vgjj((/>) explicitly depends on time t and this dependence is not surpris¬ 
ing for a physical quantity calculated in a time evolving background. The time in (19) is naturally 
specified by the correlation function of interest that is aimed to be determined in the beginning. 
Indeed, the time integrals in (14) or (18) are also limited by this initially assigned time parame¬ 
ter. Unlike the in-out computations in flat space that extends in the whole timeline, in the in-in 
framework one may imagine integrating out the “modes” from some initially prescribed time to to 
the time of interest t. 

In carrying out the x path integral above no approximation is utilized. Consequently, given the 
initial interaction potential (2), the result (19) is exact. On the other hand, one may attempt for 
an alternative but perturbative path integral calculation that can be organized in powers of the 
external cj) field, where the corresponding series can be pictured as in Fig. 1. The effective potential 
(19) equals the sum of this infinite series and conversely the infinite series can be recovered from 
(19) by expanding it in powers of (f). 

^ The computation of the effective potential can also be done by first introducing the Fourier modes in the path 
integral, which factorizes (6) for each momentum mode k. This yields the determinant of an ordinary differential 
operator for each k, which can be calculated using (A.17). The final result involves a sum over the momentum 
modes that appears as the momentum integral in (19). 
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For gcj) <C H, the (j) dependence in the mode equation (16) becomes negligible. In this limit, one 
hnds Vgjj ~ which is identical to the one-loop Hartree approximation applied to the 

original potential (2). 

It is instructive to use the formula (19) in flat spacetime. Using the mode function 


Mfc = 


\/%Jk 


^-icokt 


where 01 “^ = k? + fh? + sg^cjfl, in (19) and carrying out the s-integral, one finds 

1 

= dkk^ 


\/k‘^ + — k 


In [15], the bare Coleman-Weinberg potential for the interaction Xcffl/Al is given by 


vHV = 


1 


/ 


d^kf 


■ In 


1 + 


X(f) 


2k‘^ j 


( 20 ) 


( 21 ) 


( 22 ) 


2J (27r)2 

In obtaining (22), one dehnes a fluctuation field as 4> = (pB + d(p, which gives a quadratic potential 
and integrates over 5(j)- By comparing this quadratic potential with the potential we 
are dealing with, i.e. one sees that (22) must be identical to (21) after the replacements 

4>b ^ (p and A — )■ 2g^ (the helds that are integrated out are 5(p and y, respectively). Indeed, 
defining kE = {k,kfl} (so that /c|; = k.k + k^ = k"^ + k^) and carrying out the k^ integral in (22) 
exactly yields (21). 

The form of the bare effective potential (19) is suitable for adiabatic regularization. The adia¬ 
batic mode function is defined as 


= 


1 


V‘^^k 


and the mode equation (16) implies 


exp —i / 


snl ink 


(23) 


where 


02 _ , ,2 , ^ "‘fc 

“*-“* + 401 “2 S5’ 


9 k _ 9 9,9 

(jjt = ^ + m + sg (p . 


(24) 


(25) 


Note that (23) obeys the Wronskian condition identically. One may now solve (24) iteratively in 
the number of time derivatives starting from the 0th order solution Truncating the 

iterative solution at any desired order^ defines an approximate adiabatic vacuum of that order. 


In some cases, there can be an inherent ambiguity in defining the order in this iterative solution, see [29]. 



Since (19) is expected to be quadratically divergent (this is the case for the Bunch-Davies 
vacuum), it is enough to subtract the second order adiabatic terms from (19). A straightforward 
calculation gives the renormalized effective action as 


K//(0) = 


1 


2a{t)'' 


2 i.2 

■9 (t> 


d?k 

(27r)3 


ds 


\9k[t,s,4>\\ - 




(26) 


where 




1 

2 cufc 


9i72 3iJ 
8 4 8 4 a?ujl ' 2 a?ujl 


1 HP 1 772^2 

+ X' 


(27) 


Note that Wfc, which is defined in (25), depends on the parameter s. Eq. (26) is the main result of 
this section. 

The subtracted terms in (26) can be thought to arise from a counterterm potential 5V{(j)) 
appearing in the bare action. It can be found as 


SVicP) 


-aig 


- 2aig‘^rP(j)^ - 


0277 + a^H 


(28) 


where 

IQp Jq (P + 1 ) 1/2 ’ 

_ 9 r°° Pdk 1 P dk 5 7°° 7® dk 

(72 + 1)3/2 + (72 + 1)5/2 " ^ (72 + 1)7/2 ’ 

3 /■“ 72 dk 1 P dk 

(72 + 1)3/2 " ^ (72 + 1)5/2 • 


(29) 

(30) 


This interpretation justifies the adiabatic subtraction terms in (26) since the regularization pro¬ 
cedure can be recast as a standard renormalization method (see [30] for a similar renormalization 
treatment involving adiabatic subtractions). 

One may see that in the flat space limit (26) gives Veff{(j)) = 0. This can be viewed as a 
prescription yielding an unambiguous renormalized result. Namely, the inherently present finite 
renormalizations are fixed in such a way that as the expansion of the universe is switched off the 
quantum effective potential vanishes. Thus, the adiabatic regularization physically captures the 
quantum effects directly related to the cosmic expansion or “particle creation”. 


III. APPLICATIONS TO INFLATION 

It is straightforward to apply the above formalism to any given cosmological background and 
in this section we consider the standard single field slow-roll inflationary scenario. Let (j) be the 
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inflaton field and x be another scalar field coupled to cj) in the form (2). Such couplings are necessary 
for the decay of the inflaton during reheating and the potential (2) is a plausible alternative, see 
e.g. [2, 3]. Although, at the classical level the x field does not affect the background dynamics, one 
may wonder how the quantum x fluctuations modify the classical inflaton potential u((/>) through 
the coupling (2). For large field models like chaotic inflation, where the value of the inflaton field 
background exceeds the Planck scale, the quantum corrections might be important. 

In this computation, the slow-rolling of the inflaton field gives sub-leading corrections and 
therefore we set a = exp(Ht). The Bunch-Davies mode function of the x held in de Sitter space 
obeying (16) is given by 

s, </.] = {ke-^^/H) , (31) 


where is the standard Hankel function of first kind and 


'9 sg'^cf)^ + fh? 
- H^- 


(32) 


Using the mode function in (26) and determining the adiabatic subtraction terms for de Sitter 
space from (27), we obtain 




du / ds 






u\ 


^ 9 V? 

1 + X—+ 




2uj{u) \ 8a;(u)2 2a;(u)^ 8uj{u)^ J \ 


, (33) 


where 


/ n 2 2 , + fri^ 

UJ[u) = u + 




(34) 


Eq. (33) gives the fully renormalized exact quantum effective potential for the inflaton field that 
arises from its interaction with the x held with the potential (2). In de Sitter space, Veff{4>) turns 
out to be time independent due to the underlying symmetries of the background. 

If m 7 ^ 0, the u-integral is convergent near u = 0 for any value of (j). When rh = 0, i.e when 
X is massless, one may see that as (p ^ 0 the integrand in (33) vanishes and thus Vf,ff{(p) — )• 0. 
On the other hand, for nonzero cj) the u-integral is again convergent as u —)■ 0 since Re(i/) < 3/2. 
Moreover, the adiabatic regularization guarantees UV finiteness and thus the convergence as u —)■ 
oo. Therefore, (33) is completely safe both at IR and UV, corresponding to u —)■ 0 and u —)■ oo 
limits, respectively. 

As mentioned above, (33) is the quantum correction to the classical inflaton potential due 
to its coupling to the x held and it can be used in Einstein’s equations to modify the classical 
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0.01 

0.1 

1 

2 

4 

20 

100 

200 


0.005 

0.03 

0.06 

0.07 

0.07 

0.08 

0.07 

0.06 


TABLE 1. The numerical values of (33) for various field values in units of the Hubble parameter H. 

evolution (although the corrections turn out to be small as we will show below). This effect is 
somehow complementary to the one obtained by the so called tadpole method, where one defines 
a classical background and a fluctuation field, and the evolution of the background is fixed by the 
vanishing of the tadpole of the fluctuation field, see e.g. [31~33]. Indeed, one may see that at one 
loop the tadpole method applied to the interaction potential in (2) is equivalent to the Hartree 
approximation, where appears as an effective mass term for the inflaton. At higher loops, 

the tadpole method necessarily includes the inflaton loops which are completely absent in (33). 

In a single field model, it is natural to assume that m H, since otherwise the y field is 
expected to alter the inflationary background evolution unless the initial conditions are finely 
tuned. Neglecting m in (33), one sees that the effective potential becomes times a dimensionless 
function For gcf) ~ H, the function / takes values of order unity and its dependence 

on its argument is somehow weak. In Table I, we numerically evaluate Veff for various values of 
g(l)/H. Since 14// ^ and the background energy density is of the order of the corrections 

induced by 14 // on the cosmological observables are suppressed by the factor H'^/Mp. 

Let us now discuss a similar problem involving the cosmological perturbations. We focus on the 
curvature perturbation which can be introduced in the metric as 

ds'^ = —dt^ + a{t)'^ (35) 

As before, we treat C as an external field and integrate out the y field in the path integral to obtain 
an effective action. In the presence of the curvature perturbation (in the = 0 gauge) the metric 
(35) must include nontrivial lapse N and the shift iV* functions that depend on ( [34]. Nevertheless, 
we take C to be an (off-shell) superhorizon perturbation and apply a derivative expansion to the 
effective C action to determine an effective potential. In that case () dependent factors in both N 
and become negligible at the leading order since they always contain derivatives. We take y to 
be a minimally coupled massless field propagating in the metric (35) where the action becomes 

5 = ^ y d‘^xa^e^‘^ y^ - ^e~‘^^{dxf ■ 


(36) 
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Note that there is a shift symmetry in the action that is given by 

C —5- C + (37) 

and we demand this symmetry to be preserved in our computation below. 

The path integral over the momentum variable is nontrivial since it involves the external ( 
field. After applying the following canonical transformation 

exp [3C/2] y, 

PfM = exp [-3C/2] P^, (38) 

which preserves the measure 


DxDP^ = D^DP^, 


(39) 


the path integral over becomes independent of C and decouples. As a result, one is left with an 
integral over ^ with the action 


5 = 


/ 


d'^x 


-2C 


■ 2 I / 3 Tj I 9 tt 2 \ 2 ^ 

+( 2 ^+ 4 ^ 


iduf 


(40) 


where the derivatives of C, are again neglected for consistency with the original action (36). This 
final path integral yields the functional determinant of the operator in (40). It can be seen that 
the in-in calculation is again separated into two identical + and — branches, as in the previous 
section. 

Even for a slowly changing external held it is very difficult to calculate this functional de¬ 
terminant exactly. On the other hand, the constant C conhgurations are pure gauge due to the 
shift symmetry and the exactly calculable determinant does not give any information about the 
effective action. Therefore, we focus on the curvature perturbations obeying 


C«i7, 

diC < /c*C 


(41) 


for some A:* so that the path integral can partially be carried out as follows: One may see that the 
time derivatives of C, always appear with H in the combination (^ -|- H), thus the hrst condition in 
(41) allows one to neglect C, in the action. Decomposing the path integral measure in the momentum 
space as 




(42) 




k k>k* 
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and integrating out the modes with k > the spatial derivatives of C can also be neglected in 
the action. As a result, ^ can effectively be treated as a constant in such a partial path integral of 
“high energy modes”. For an actual superhorizon curvature perturbation both conditions in 
(41) are satisfied once /c* is taken to be the comoving horizon scale determined by the metric (35). 
Therefore, by setting 

/c* = (43) 


it is possible to obtain an effective potential for (sh by integrating out the subhorizon y modes. 

In normalizing this path integral, one must preserve the original shift symmetry (37) (this is 
necessary due to the way ( is introduced in the metric (35)). For that it is convenient to calculate 
the derivative of the functional determinant from (A. 12) (with the s derivative replaced by a 
(sh derivative), which yields the derivative of the effective potential. The undetermined additive 
constant of integration can uniquely be fixed by demanding the symmetry (37). A straightforward 
calculation using the formulas derived in the previous section gives an effective bare potential for 


Csh as 


r5/(a)' = / 

a Jk>kt 


; 2 I 12 

(27r)3 ^ ’ 


(44) 


where the prime denotes derivative with respect to argument and the mode functions obey 


^2 

/ffc + -^k'k- 

k^kfj*k - dAfj'k = i- 


9 


3 • 


+ -H 


l-^k — 0 , 


(45) 


We note that (44) is exact to the leading order in the derivative expansion. 

In the Bunch-Davies vacuum, the bare potential (44) has a quartic UV divergence that can be 
cured by adiabatic regularization by subtracting the terms up to fourth order. This gives 


1 


VeffiCshY = -^e 


lk>k, 


(fk 

(27r)3 




w- 


(4) 


(46) 


where 


(4) 

k'ad 


— ^Csh 


a 


1 + + oo) + (^3d"^ + 3add — 5a^da^^^ — 


p4C., 


(47) 


One may see that in the exact de Sitter space (46) vanishes,^ which is a standard result in adiabatic 
regularization (see e.g. [22]). Therefore, the slow rolling of the inflaton must be taken into account 

^ In applying the adiabatic regularization to massless fields, one must actually start with a massive field and then take 
a limit. In this procedure, one discovers some terms that survive “mass going to zero” limit when the momentum 
integral is unconstrained and this is how the standard trace anomaly is recovered in the adiabatic regularization 
[35]. In our case, since the momentum integral in (46) has already an IR cutoff, this procedure gives no extra term 
in the massless limit. 
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e 

10"! 

10-1 

10-1 

10-2 

10-2 

10-2 

10-3 

10-3 

10-3 

a 

1 

10 

100 

1 

10 

100 

1 

10 

100 

c 

10-3 

4 X 10-1 

10-4 

1.1 X 10-4 

1.0 X 10-4 

0.9 X 10-4 

1.10 X 10-3 

1.08 X 10-3 

1.07 X 10-3 


TABLE II. The numerical values of the coefficient C given in (52) for different values of e evaluated at 
different times labeled by the scale factor a. Note that for this set C ~ 10“^e and the pre-factor is related 
to the overall number l/(67r^) appearing in (52). 


and here we take the background to be the space with constant deceleration obeying 

H = -eH^, (48) 

where e is a constant that is not necessarily small. Introducing the conformal time adr] = dt, the 
scale factor can be expressed as 

a = (49) 


where the constant Hq can be identified with the inflationary Hubble scale when e <C 1. The 
Bunch-Davies mode function obeying (45) is given by 



and (46) becomes 


3-e 

“ 2(1-e) 


(50) 


Veff{Csh) = CH^e^i^^ 


(51) 


where 


C = - 


1 

Gvr^ 


du 


IT 




2u 


2-e 3e(2-e)(3-2e) 

2(1 — 8(1 — 


(52) 


In fixing the integration constant in (51) we demand the action to have the symmetry (37). Note 
that the coefficient C is time dependent, but this dependence is mild in the slow-roll regime. 

One may see that C —)■ 0 as e —)• 0. Thus, C ~ e for e <C 1 (see Table II) and in the slow-roll 
approximation (46) becomes 


Veff{Csh)-eHle^^^K 


(53) 


It is important to emphasize that (53) can only be used for superhorizon modes and the corrections 
induced by (53) occur after horizon crossing. Physically, one may imagine them to arise due to 
the scatterings of the super horizon Qgh with the subhorizon (virtual) x fluctuations. The theorems 
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proved in [36, 37] regarding the conservation of superhorizon ( in quantum theory are not applicable 
here due to presence of the extra scalar field y. 

One may calculate the tree level corrections to the power spectrum and non-gaussianity that 
follow from (53). The quadratic piece coming from the expansion of the exponential yields a mass 
term that modifies the power spectrum and the cubic term yields a nonzero Jnl- The relative 
magnitude of the correction to the power spectrum and /nl are given by 


N 


M2 


(54) 


where N is the number of e-folds from the horizon crossing (of the mode of interest) till the end 
of inflation (with the start of reheating loop corrections may still be important as discussed in 
[38-40] but the whole calculation takes a different form). This slight enhancement is similar to 
ln(a) factors appearing in loop corrections, see e.g. [16]. 

Finally, it is also possible to work out an effective potential for the inflaton that arises from inte¬ 
grating out its own quantum fluctuations. To see that, one may introduce the standard generating 
functional® 


omj] - 


= / D(^e 


iScl[4>]+i S J<P 


(55) 


and the quantum effective action 

T[cl>d] = W[J]- J Jckci. 

Using ST/5(pci = — J, (55) can be written as 

which is an exact expression for T. Expanding the quantum effective action as 

T = S,i + S^^\ 

where denotes a first order (one-loop) correction, defining a new integration variable x = 
and keeping only the leading order terms one finds 


(56) 


(57) 




6^ Sr 


I I U-Jcl 
O / X .. .. X 


~ cl' 


(58) 




(59) 


It is straightforward to repeat the analysis below for the in-in case. 
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This gives the one-loop effective action in terms of the functional determinant of the (second 
order) fluctuation operator. Starting with a classical inflaton potential v{(j)), the adiabatically 
regularized one loop quantum potential is given by (26) with the replacement 

9^ ^ (60) 

where the prime denotes derivative with respect to argument. A similar analysis can be done for 
the curvature perturbation 


IV. CONCLUSIONS 

In this paper, we use mode functions to express the functional determinants that give one-loop 
effective potentials of scalar fields propagating in cosmological backgrounds and apply adiabatic 
regularization to cure the resulting infinities. It turns out that in this problem the adiabatic 
regularization can be cast into a well defined renormalization procedure corresponding to the 
specific counterterms added to the action. The finite renormalizations are automatically fixed in 
such a way that the effective potential vanishes when the expansion of the universe is turned off, 
thus it captures the quantum effects related to the expansion. 

The same method can also be applied to derive an effective potential for the superhorizon 
curvature perturbation (^sh encoding its interactions with virtual subhorizon modes circulating 
in the loops. For a given mode, the effective potential can be used starting from the horizon 
crossing time till the beginning of reheating. In that case, it is difficult to interpret the adiabatic 
regularization by counterterms added to the action since the momentum integral is cut off at the 
horizon crossing time (see e.g. slow-roll parameter dependent terms in (52)). In any case, the 
loop effects involving cosmological perturbations already contain non-renormalizable gravitational 
interactions and consequently the whole problem should be viewed in an effective field theory 
approach. 

The corrections induced by the effective inflaton potential on the background evolution and 
effective Csh potential on the cosmological observables like power spectrum and non-gaussianity 
turn out to be small even in the large field models like chaotic. These findings support the general 
view that inflation can fully be understood in a semiclassical approach and quantum backreaction 
effects can safely be neglected. 
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Appendix: The functional determinant of a time dependent in-in operator 


We consider the Hilbert space of doublets of time dependent functions 

m = 


(P+it) 

(p-{t) 


(A.l) 


which are defined in the interval (tajtb) with the inner product 

<$i|$2>= / dt^J{t)^2{t) = [ dt [(pf {t) + (t)] . (A.2) 

Jta Jta 

We assume that these functions obey suitable boundary conditions so that the one parameter 
family of operators L^, which have the following form 


L 0 
0 -L 


+ s 


(A.3) 


H+(f) 0 

0 

are self-adjoint with discrete spectrum. Here, L denotes a generic second order (ordinary) differ¬ 
ential operator, s is a real parameter and are smooth functions. Our aim is to calculate the 
ratios of the two determinants in this family, i.e. detLi/detLo- 
For a given parameter s, let 


Vnis) = 


ynit,s) 

to denote the eigenfunctions of the self-adjoint operator (A.3), i.e. 


(A.4) 


Ls2/n('S) — ^n(s)?/n('S). 


The eigenfunctions can be taken to be orthonormal 


<yn{s)\ym{s)>= 5n 


(A.5) 


(A.6) 


and they also form a complete set of basis functions in the Hilbert space satisfying the following 
completeness relation 


'^ynis,t)y'^is,t') 


5{t — t') 0 

0 6{t — T) 

The Green function with the prescribed boundary conditions is defined by 


(A.7) 


L,G, = 


d{t — t') 0 

0 6{t — t') 


(A.8) 
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We assume that has no zero modes and the Green function is uniquely defined. The matrix 


entries, which can be introduced as 


Gsit,t') = 


are fixed by the eigenfunctions yn{s) as 


G++{t,t') G+-{t,t') 
Gj+{t,t') Gj-{t,t') 


G.(t,0 = E 


yn{s,t)y^{s,t' 

^n(s) 


(A.IO) 


As discussed in [26], one may use this identity in constructing standard in-in propagators in flat 


space. 


Formally, the determinant can be written as 


det Ls = exp [Tr In L^j , 


(A.ll) 


where the trace is defined in the Hilbert space (A.2). Differentiating the above expression with 


respect to s and noting that ^ = G* one finds 


d , , ^ ^ dLg 

— lndetLs = Tr G^—— , 
as as 


(A.12) 


which can be used to express the determinant an operator in terms of its Green function (see e.g. 
[4-9]). 

To get a more rigorous derivation of (A.12), one may follow a simple method presented in [14]. 
The determinant of the operator, which is defined as 


det Ayi(s), 


(A.13) 


generically diverges. However, the ratio of the two determinants like detLi/ detLo is usually well 
defined and finite. Differentiating the eigenvalue equation 


K(s) =<yn(s)lLslyn(s)> 


(A.14) 


with respect to s (an s-derivative is denoted by a prime), one may find 

0 

A„(s) =<|/„(s)l lyn(s)> . 

0 -Q- 


Then, 


(A.15) 


^ M' \^f*'’jJyn(i,s)y+(t,s) ^ ynit,s)y^it,s) 

Eln[A.(s)] dt ^ n it) -(f)J, 


(A.16) 
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where the components of the eigenfunctions are introduced in (A.4). Integrating this last equation 
with respect ot s and using (A.9) and (A. 10), one finally obtains 
det Li 


In 


= [ ds [, (A.17) 

Jo Jta 


det Lq 

which expresses the ratios of the determinants in terms of the Green function. 

One nice feature of the construction of [14] is that it can be readily generalized to partial dif¬ 
ferential operators, unlike the case with the Gel’fand-Yaglom theorem (see e.g. [41, 42]). Consider 
the doublets of functions 

0+(x^) 

4>~{x^) 

which are defined in an n-dimensional space parametrized with the coordinates and endowed 
with the following inner product 


4>(x) = 


(A.18) 


<4>i|4>2>= J (x)4>2(x). 


(A.19) 


Consider a self-adjoint operator of the form (A.3) where L now denotes a partial differential operator 
and n^(t) —)■ x^). It is easy to repeat the above steps to show that 

det Li 


In 


= J ds J dP'x [G^+(x; x)n’''(x) — (x;x)ll (x)] 


(A.20) 


L5Gs(x;x') = 


(A.21) 


det Lq 

where the Green function is defined as 

(5"'(x —x') 0 

0 5'^{x — x') 

In the cosmological setting of our interest, one has x^ = (t,x*) and L is fixed by the covariant 
Laplacian as in (10). In that case, the spectrum of the operator is usually continuous. However, 
one can initially take the spatial coordinates x* to be periodic with period a giving a discrete set of 
momenta 2'nrd ja with integers rd. One may then take a —)■ oo limit that would replace the discrete 
Fourier modes with the continuous Fourier transform, which would justify the use of (A.20). 
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